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We study numerically a stochastic differential equation describing an interface driven along 
the hard direction of an anisotropic random medium. The interface is subject to a homogeneous 
driving force, random pinning forces and the surface tension. In addition, a nonlinear term due 
to the anisotropy of the medium is included. The critical exponents characterizing the depinning 
transition are determined numerically for a one-dimensional interface. The results are the same, 
within errors, as those of the "Directed Percolation Depinning" (DPD) model. We therefore expect 
that the critical exponents of the stochastic differential equation are exactly given by the exponents 
obtained by a mapping of the DPD model to directed percolation. We find that a moving interface 
near the depinning transition is not self-affine and shows a behavior similar to the DPD model. 
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I. INTRODUCTION 

The driven viscous motion of an clastic interface in a 
medium with random pinning forces is relevant for the 
understanding of various problems in condensed matter 
physics Examples include the ordering dynamics of 
an impure Ising magnet after a quench below the critical 
temperature [^, wetting immiscible displacement of one 
fluid by another in a porous medium H, and pinning 
of flux lines in typc-II-supcrconductors^,^. In recent 
years, studies of fluid imbibition in paper and of 
flameless paper burning j|] have been carried out, and 
observable interfaces allow a direct comparison with the- 
oretical predictions. 

Common to all of these problems is a competition be- 
tween smoothening due to the surface tension and rough- 
ening due to the interaction with the random pinning 
forces of the medium. Further, there is a competition 
between the driving force and the pinning forces, result- 
ing in a depinning transition. 

On a coarse-grained level, it is expected that the dy- 
namics of the interface can be described by the following 
continuum equation of motion. 



Vn{r, t) = -iK{r) + F + T]{f) + n ■ VV{r). 



(1) 



Here, w„(r, t) is the normal velocity of the interface at 
position r. The surface tension generates a term pro- 
portional to the total curvature K(r) = — V • n, where 
ft is the normal vector on the interface at r. The coef- 
ficient 7 measures the stiffness of the interface. F is a 
homogeneous driving force. The last two terms, 77(7^ and 
n ■ W{f) represent random- field and random-bond dis- 
order, respectively. The random forces r](f) and the ran- 
dom potential V^f) are short-range correlated in space. 



Equation (Q) is considerably simplified by restricting 
ourselves to an almost planar interface without over- 
hangs. A coordinate system r = (x, h) can be introduced, 
so that the interface position is given by a single-valued 
function h(x,t). The dimension of x is denoted by d. 
Equation (|l|) becomes j|] 



1 dh{x,t) 



7 



V9 dt 



9 



1 



VhVV{x,h) - 



dV{x, h) 



dh 



(2) 



where g = 1 + (V/i)^. 

For sufficiently large values of the driving force F, the 
interface grows continuously. However, for smaller values 
of F, growth on some regions of the interface can come to 
a halt, due to the interaction with the quenched disorder. 
We say that these regions of the interface have become 
pinned. As the rest of the interface continues to grow, 
the pinned regions can be dragged over the pinning bar- 
riers by neighboring moving regions. Then, the formerly 
pinned regions advance quickly, which can be considered 
as an avalanche ]To|-p^. 

The maximum linear size, ^, of the pinned regions di- 
verges when F approaches its critical value F^, 



\F-Fr 



(3) 



The threshold Fc is the critical point of a dynamical 
phase transition, and ^ the corresponding correlation 
length. The role of the order parameter is played by the 
mean velocity, v = lim^^oo.L-^oo dh/dt. {L is the system 
size and the overbar denotes the spatial average over x.) 
The velocity is zero for F < Fc, and increases as 
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(4) 



for F > Fc- On length scales / 3> C pinning can be 
neglected and we can therefore replace the argument 
h in the disorder terms 77(x, /i) and V(j>i,h) by vt, i.e., 
the quenched disorder crosses over to thermal noise @ . 
Then, the interface is gove rned by the Kardar-Parisi- 
Zhang (KPZ) equation In this paper, we are in- 

terested in the critical behavior on length scales I <C 
especially when ^ — > oo at the depinning transition. 

The global interface width = ( h{x, t) — h{t) ) is an- 



other characteristic quantity of the interface. Here and 
elsewhere, (•) denotes an average over the disorder distri- 
bution. Choosing a flat interface as the initial condition, 
/i(x, t = 0) = 0,w^ scales as MK% 



(5) 



for a sufficiently large system size, L > ^. If L < ^, 
the correlation length ^ in Eq. (^) has to be replaced by 
L. ^+{y) and *_(y) are scaling functions for F > F^ 
and F < Fc, respectively. Both functions scale as 
^±(y) ^ y'^'^ for 2/ <C 1, where f3 — a/z. It follows 



(6) 



For F > Fc, pinning is irrelevant on length scales I ^ ^, 
so we can neglect pinning also on time scales t ^ . 
Thus, ^ y2/3m Jqj. y ^ where f3rn is the growth 

exponent of an interface subject to thermal noise. Below 
threshold, the interface becomes pinned and ^'-(y) — 
const. Using the scaling of 'i'±{y) for y <C 1 and Eq. (|^) 
we obtain 



rt2/3™ {t:$>e,F>Fc) 
\e" {t:^e,F<Fc). 



(7) 



It has been shown that the critical exponents fulfill an 
exact scaling relation |l4|, 



^{z — a). 



(8) 



II. 



THE MODEL FOR ANISOTROPIC 
DEPINNING 



To further simplify the equation of motion (|), we as- 
sume that the typical gradients Vh are small on large 
length scales, so that the roughness exponent a is smaller 
than one. This assumption has to be compared with the 
final results. When expanding ~ 1 — (V/i)^/2, non- 

linear terms proportional to (V/i)^ are generated in Eq. 
(^. A natural question is whether these terms are rel- 
evant at the depinning transition. A term (V/i)^ with 
a positive coefficient on the right hand side of Eq. (||) 
would give a nonzero contribution to the driving force 
for any rough interface. This contribution increases when 



imposing a global tilt of the interface. Thus, the thresh- 
old Fc becomes a function of the average orientation of 
the interface. This is reasonable for anisotropic systems 
but not for isotropic ones. For interfaces in an isotropic 
environment, it can indeed be shown that the nonlinear 
terms generated by expanding 1 / in Eq. (||) are irrel- 
evant close to the depinning transition pT[ |. If however, 
the medium is anisotropic, a term of the form A(V/i)^ 
can be relevant even for the case f — > 0. In fact. Tang, 
Kardar, and Dhar |^ argued that if W and dV/dh are 
differently distributed, the term A(V/i)^ is generated un- 
der coarse graining. More generally, if the system is 
anisotropic in the sense that the threshold Fc depends 
on the average orientation of the interface, A(V/i)^ is the 
only relevant term that can change the universality class 
of the depinning transition Motivated by these ob- 
servations, we consider the following equation of motion. 



dhjx, t) 
dt 



-fV^h + A(V/i)2 +F + r]{x, h). 



(9) 



For A > the threshold Fc has a maximum for an inter- 
face without tilt, i.e^^an interface with periodic boundary 
conditions in Eq. is driven along the hard direction 
of the anisotropic medium . 

Equation (W is the model we study in this paper. 
For simplicity, we restrict ourselves to interface dimen- 
sion d = 1 and consider only random-field disorder 
ri(x,h). It was shown by Narayan and Fisher [ pl| for 
isotropic systems that random fields and random bonds 
give rise to the same critical behavior. This has been 
supported by numerical simulations of interfaces subject 
to random-bond disorder [ p^|Jl9|] . The random forces are 
assumed to have zero mean and short-range correlations, 
{■q{x',h')r]{x' + x,h' + h)) = S{x)A{h), where A{h) de- 
creases exponentially for \h\ greater than a microscopic 
cutoff. 



Previous Results 

The case A = in Eq. (||) was first investigated by 
Feigel'man |20|. Significant progress has been made by 
a functional renormalization-group treatment |pl] , p^ , pl| 
and by extensive numerical simulations p3| , pl} -p5|. 

The results for the anisotropic case A > are less well 
established. It was first suggested in Ref. ||2^] that Eq. 

is in the same universality class as the "Directed Per- 
colation Depinning" (DPD) model |,|l|l. For d = 1, 
directed percolating paths of pinning sites stop the inter- 
face. Thus, the roughness of the pinned interface is given 
by the scaling of the directed paths and the correspond- 
ing roughness exponent is a ~ 0.633 §,0. The other 
critical exponents of the DPD model can also be obtained 
by a mapping to directed percolation; in d = 1 the dy- 
namical exponent z — 1, the correlation length exponent 
ly ~ 1.733, and the velocity exponent 9 ~ 0.636. 

Amaral, Barabasi, and Stanley |^ measured the tilt 
dependence of the velocity for several versions of the DPD 
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model and found a behavior consistent with the relevance 
of the term A(Vfe)^ at the depinning transition. Galluccio 
and Zhang Q simulated a self-organized version [ p9|j30| 
of Eq. (H) and obtained a ~ 0.63, thereby supporting the 
conjecture of Rcf. Recently, Olami, Procaccia, and 
Zeitak |25| argued that the slopes of the pinned surfaces 
in Eq. are bounded and that therefore Eq. could 
belong to the directed percolation universality class. 

However, Csahok et al. performed the only direct 
simulation of the continuum equation (^) and came to 
a different conclusion: Their numerical values for the 
exponents a and /3 in d = 1 are in agreement with a 
scaling theory for Eq. w hich yields a = (4 — (i)/4 
and /3 = (4 — d) /JA + d) |31| . Yet another proposal was 
made by Parisi 1 32 1 . He argued that /3 = (4 — (i)/4 and 
supported this value in c? = 1 by a simulation of a lat- 
tice model, which was assumed to be in the universality 
class of Eq. (^. Problematic is also the interpretation of 
Stepanow's renormalization-group calculation in d = 4— e 
dimensions js^ . The extrapolation to d = 1 gives the re- 
sults a ~ 0.86, z~l,iy~ 1.2, and 9 ~ 0.16 Is^]. 



Aim of the Paper 

To resolve these discrepancies we carefully determine 
in this paper the critical exponents by large-scale simu- 
lations of Eq. (H) for interface dimension d = 1. To this 
end we carry out a numerical integration of the equation 
of motion with a continuous height variable as well as 
a simulation of an automaton model where the height 
variable takes integer values only. In addition to the pre- 
viously measured exponents a and /3 we also determine 
the exponents v and 9 to strengthen our conclusions. We 
find that the numerical values for all critical exponents 
are in agreement with the suggestion that Eq. (^) and 
the DPD model belong to the same universality class. 

In addition, we investigate the interface roughness for 
F > Fc- It is shown that a moving interface is not self- 
afhne and that the behavior of the roughness is very sim- 
ilar to that of the DPD model. 



III. NUMERICAL METHODS 

A. Continuum equation 

First we simulate Eq. (^) with a discretization of 
the transverse coordinate only, cc — > i, h(x,t) hi(t) 
(with 1 < i < L). The random forces r]i{h) are cho- 
sen as follows: Each integer position hi on a square lat- 
tice is assigned a random number 77 between zero and 
one. For non- integer hi the forces r]i(h) are obtained by 
linear interpolation |Q. Finally, the /i-coordinates of 
r]i{h) in each column i are shifted by a random amount 
< Si < 1, i.e., rii{h) rii{h + si) 



At t = the interface is flat, hiit = Q) = Q. The inter- 
face configuration aX t -\- At is calculated simultaneously 
for all i using the method of finite differences j35| , 

h,{t -t- At) = h,{t) + At{7[/i,+i(t) + h,^i{t) - 2h,{t)\ + 

+ X[h,+i{t) - /i,_i(t)]' + gm{hi)]. (10) 

Periodic boundary conditions are used and 5 is a parame- 
ter measuring the strength of the disorder. We choose the 
parameters 7 = 5, A = 1, 17 = 3 |^ and use At — 0.04, 
for which Eq. (RG) is found to be stable. We checked 
that simulations with At = 0.01 yield consistent results. 

B. Automaton model 

Since the simulations of the continuum equation (|l0| ) 
are conrputationally expensive, we also study a lattice 
model of probabilistic cellular automata, which al- 
lows to determine the critical exponents more effectively. 

The automaton model is defined on a square lattice 
where each cell [i, h] (with 1 < j < L) is assigned a ran- 
dom force rji^h which takes the value 1 with probability p 
and rji^h = —1 with probability \ — p. During the motion 
at a given time t the local force 

f,{t) - 7[/i,+i(t) + h,_^{t) - 2K{t)\ + 

+ \[hi+i{t) - hi-i{t)\^ + grii^hi (11) 

is determined for all i. The interface configuration is then 
updated simultaneously for all i 122]: 



h,{t + \) = h,{t) + 1 
h^{t + 1) = h,{t) 



if /. >0 
otherwise. 



(12) 



The difference p—{l~p) ~ 2p—l determines the driv- 
ing force. For simplicity, we use the density p as the tun- 
able parameter, and the depinning threshold is denoted 
by pc- The results shown in this paper were performed 
with the following parameters, 7 = 10, A = 1, and g = 20 
|6|. 

The growth rule specified by Eqs. (11) and ( [T^ ) can 
be derived from the continuum equation ( 9|) by temporal 
and spatial discretizations. In addition, a simple two- 
state random force r]i^h — ±1 is used. The discretization 
implies that the critical slowing down close to the thresh- 
old is reduced. This can be seen as follows. For w — > 0, 
the local force for most of the interface elements hi (t) in 
Eq. (|l^) is almost zero. Nonetheless, in the finite differ- 
ence approximation Eq. (p^), all hi{t) and rji^^i have to 
be updated at each time step. On the other hand, only a 
subset of values hi{t) and rji^hi are updated at each time 
step in the cellular automata model. 
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IV. NUMERICAL RESULTS 



A. Roughness at threshold 



In this section we determine the critical exponents a 
and P, defined in Eqs. (^) and (^), respectively. First 
we find the threshold value Fc and Pc from a measure- 
ment of the interface width ■w'^{t) for different values of 
the driving force. Since the correlation length increases 
when F ^ Fc, the range of the scaling regime w'^{t) ~ t'^^ 
also increases (see Eq. (^). The threshold is estimated 
as the value where the power-law scaling holds for the 
longest time interval. This method allows to determine 
the threshold very accurately. 
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FIG. 1. (a) Scaling of the interface width w^it) 
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at threshold. Interfaces of size L = 16384 at — 0.1511 
were simulated for the continuum equation (10), and at 
Pc ~ 0.6631 for the automaton (L = 262144). The data 
were averaged over 40 independent disorder distributions. 
The statistical uncertainties are smaller than the size of 
the symbols. The asymptotic slopes yield (3 — 0.64 ± 0.02 
for both lines, (b) The corresponding effective exponents 
f3{t) = log[w'^{t)/w^{t/2)]/log4. The value /? ~ 0.633 of the 
DPD model is shown as the dashed horizontal line. 



In Fig. |l|a the width w'^{t) at the estimated thresh- 
old is shown in the transient regime, t <C 5- In 
Fig. Jl|b the corresponding effective exponents /3(t) = 
log[w^t)/ii;^(t/2)]/log4 are shown. Plotting the effec- 
tive exponents P{t) sensitively shows the quality of the 
power-law scaling. From Fig. |l| we conclude that 
f3 = 0.64 ± 0.02 for both the continuum equation ( |lO| ) 
and the automaton (Eqs. (|ll|), ([l^)). This is in very 
good agreement with the value (3 ~ 0.633, obtained by 
mapping the DPD model to directed percolation . 

A convenient estimate of the roughness exponent a 
can be obtained by measuring the equa l-time correlation 



function C (r, t) 



faces. In Fig 
and F < Fr is s 



[hi+r{t) — hi{t)]'^) of pinned inter- 



a, the scaling C^(r) ^ r^" for r <C ^ 
lown. The corresponding effective ex- 
ponents a(r) — log[C^(r)/C^(r/2)]/log4 are plotted in 
Fig. and Equation (|o|) and the automaton yield 
the same result, a = 0.63±0.01, again in agreement with 
the prediction of directed percolation, a ~ 0.633. Using 
the numerical result for (3, we see that the simulations 
are also consistent with z = a/ (3 = 1 P,p^. 



B. Roughness above threshold 



In this section we consider the steady-state behavior 
of the roughness for moving interfaces at driving forces 
F > Fc such that ^ < L. We are interested in the large- 
time limit, t ^ when the instantaneous velocity of 
the interface fluctuates around its mean value v. 

Recently, it was proposed that the roughness of mov- 
ing interfaces in the DPD model exhibits scaling with 
exponents different from the critical ones |Q. This 
is in disagreement with Ref. |jl^, where it was argued 
that on scales I < ^, the moving interfaces are not 
self-affine, because they consist of pinned regions with 
a ~ 0.63 and laterally moving regions with roughly lin- 
ear slopes (a ~ 1). As a consequence, different mo- 
ments of the equal-time correlation function, C'(r, i) = 
{\hi+r{t) — hi{t)\i), yield different effective roughness ex- 
ponents, aq{r) = log[C"'(r)/C«(r/2)]/[glog2]. 

We measure C (r, t) for the continuum equation ( |l0| ) 
and for the automaton and find results very similar to 
the DPD model (see Fig. |b and |c). The effective 
exponents aq{r) increase with q. The reason is that the 
moving regions (large slopes) have an increasing weight 
with increasing q p7[ |. 

Another possibility to investigate the scaling of moving 
interfaces is a measure ment of the height-h eight correla- 
tion function, c'^{t) = {[hi{t + t) - hi{t)]'^), for t » 
and < L. For r ^ the height-height correlation 
function scales as (?{t) ~ r^'^, provided the interface is 
self-affine. As an illustration we consider the case A = 0. 
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Height-height correlation function for the case A = 



The growth exponent /? has been determined in Refs. 
p2|,p^ by simulations of the automaton model Eqs. (H]), 
( p^ ) with A = 0. The scaling of the interface width at 
threshold, Eq. (|), yields /3 ~ 0.88 || and /3 ~ 0.85 
p^ . Here, we measure the effective exponents in the 
steady-state regime, /3(t) = log[c^(r)/c^(T/2)]/ log4. 
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FIG. 2. (a) Equal-time correlation function C (r) ~ r " 
for pinned interfaces slightly below threshold: F = 0.149 
(L = 65536) for Eq. (10) and p = 0.6623 (L = 131072) for 
the automaton. We averaged the data over 12 independent 
runs. Both lines give a = 0.63 ± 0.01. (b) Critical effective 
exponents a{r) = log[C'^(r)/C^(r/2)]/ log4 corresponding to 
(a) for Eq. (10) (full symbols). The dashed line indicates the 
exponent a ~ 0.633 of the DPD model. In addition, the effec- 
tive exponents Uq (r) for different moments q of the equal-time 
correlation function C (r) for moving interfaces above thresh- 
old {F = 0.157, L = 32768) are shown: q = 1 (open circles), 
g = 2 (open squares), and q = 4 (open diamonds). The data 
were averaged over 30 disorder distributions so that the sta- 
tistical uncertainties are of the size of the symbols or smaller, 
(c) The same as (b) but for the automaton. The open symbols 
are the effective exponents for moving interfaces at p = 0.665 
(30 disorder distributions with L — 32768). 



FIG. 3. (a) Effective exponents /3(r) (full symbols) 
for the automaton with X = 0, "f — 1, g = 1, 
p = 0.801 > Pc ^ 0.8004, and L = 32768. The data were 
averaged over 30 disorder distributions. The dashed line indi- 
cates the critical exponent /3 ~ 0.88 [22] from a measurement 
of the global width at threshold. The effective exponents 
Paoij) (open symbols) seem to approach this value for large 
r. (b) The same as (a) but for the DPD model of Ref [17] 
with p = 0.467 > pc~ 0.461 (30 runs with L = 32768). 

From Fig. |^ we see that /3(r) is approximately con- 
stant over three orders of magnitude of r, with /3 = 
0.88 ± 0.01. The effective exponent (3oo{t) of the infi- 
nite moment of the height-height correlation function, 
Coo(t) = {inaxi{hi{t+T) — hi{t)}) ~ t^^ , is also shown in 
Fig. |k. The effective exponent Poo{t) increases slightly 
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for r > 16 and approaches a value (3oo which is consis- 
tent with the values of Refs. and obtained by 
the width at threshold. We conclude that the interface 
in Eq. (||) with A = is self-afhne with the same growth 
exponent /3 for both F = and F > F^. 



Height-height correlation function for the anisotropic case 

In contrast, the same measurement for the DPD model 
(see Fig. |b) gives decreasing effective exponents /3(t) 
and /3oo(''"), which shows that there is no scaling with the 
critical growth exponent (3 ~ 0.63. For small times t, 
the effective exponents /9(t) and I3oo{t) are roughly one, 
which can be understood from the lateral motion of the 
parts of the interface with linear slope. 
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Coq{t) seems to show a scaling for large r with j3oo — 0.63. 
It is not clear whether or not these plateaus in Fig. ^ can 
be considered as some (multi-) scaling regime. We can 
only conclude that the behavior is different from that of 
a self-affine interface which is characterized by the same 
critical growth exponent (3 for F — Fc and F > Fc a,s 
long as time scales r <C are considered. 



C. Scaling of the correlation length and the velocity 

We now proceed to determine the correlation length 
exponent v and the velocity exponent 0, which charac- 
terize the behavior of the interface close to the depinning 
transition. 




FIG. 4. (a) Effective exponents (3{t) and Poc{t) for the 
continuum equation [F — 0.157 > Fc — 0.1511, 30 runs with 
L = 32768). (b) The same as (a) but for the automaton 
(p = 0.665 > pc~ 0.6631, 30 runs with L = 32768). 

The behavior of the automaton with A > is very sim- 
ilar (see Fig. ^b). The effective exponent /?(t) for the 
continuum equation (Fig. ^) also decreases with t but 
shows a plateau at a value f3 ~ 0.82. The infinite moment 
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t \p-pf 

FIG. 5. (a) Scaling plot according to Eq. (5) for the con- 
tinuum equation. The same plotting symbol is used for data 
at a given driving force F. Above threshold, F is in the range 
0.163 < F < 1 and below threshold 0.03 <F< 0.12. For each 
value of F we simulated 50 independent runs with L — 16384. 
The best data collapse is achieved with a ~ 0.63, z ~ 1, and 
1/ ~ 1.7. (b) Scaling plot according to Eq. (5) for the automa- 
ton. Above threshold, p is in the range 0.664 < p < 0.695 (20 
independent runs) and below threshold 0.567 < p < 0.66 (100 
independent runs). The exponents for the best data collapse 
are a ~ 0.63, z ~ 1, and f ~ 1.72. 
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In Fig. scaling plots according to Eqs. (||) and (||) 
are shown. Since we already determined the threshold 
and the critical exponents a and z = a//?, we can tune 
the correlation length exponent v to achieve the best data 
collapse. The result for Eq. ( p^ ) is = 1.7 ± 0.1 and for 
the automaton v — 1.72 ± 0.03. The corrections to the 
scaling ~ C^" due to finite size effects are much 

larger for the continuum equation. Therefore simulations 
very close to the threshold are not shown in Fig. ||a and 
the error bar on the result for v is rather large. The re- 
sults for the correlation length exponent are consistent 
with the prediction of directed percolation, v ~ 1.733. 



= 0.32 ± 0.02 for the automaton. These values for 
13m are consistent with the exponent Pm = 1/3 of the 
KPZ equation ||l^. This supports the picture that the 
quenched disorder 77(2:, h) crosses over to a thermal noise 
ri{x,vt) on time scales (see Fig. ||). 

The scaling Eq. (jj) of the steady-state velocity v 
and the corresponding effective exponents 9 are shown 
in Fig. p. The result for the continuum equation is 
e ^ 0.64 ± 0.05, and 9 = 0.63 ± 0.02 for the automaton, 
again in agreement with the value of the DPD model, 
9 ~ 0.636 0. 
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FIG. 6. (a) The interface velocity 1; as a function of the 
driving force. The velocity is averaged over sufficiently long 
time intervals, so that the statistical uncertainties are smaller 
than the size of the symbols. For the continuum equation we 
obtain — 0.64 ± 0.05 (L — 65536) and for the automaton 
6 = 0.63 ± 0.02 (L = 131072). (b) The effective exponents 
eiF - F,) = \og[v{F)/v{F')]/\og[{F - F,)/{F' - F,)]. The 
dashed line indicates the critical value 9 ~ 0.636 of the DPD 
model. 



The growth exponent defined in Eq. (Q) for driv- 
ing forces F > Fc and time scales t 3> is found to 
be (3rn = 0.32 ± 0.03 for the continuum equation and 



V. SUMMARY AND CONCLUSIONS 

The results for the roughness of Eq. ^ with and with- 
out the term A(V/i)^ are summarized in Table I. 





Eq. (9), A = 


Eq. (9), A > 


F^FcV^O 
F > Fc, V = const., T < 
F>F„<»e 


(3 ~ 0.88 |2^] 

(3 ~ 0.88 
13m = 1/4 [pl 


a~/3~0.63 
not self-affine 

f3m - 1/3 



Table I. Behavior of the roughness for Eq. in d = 1 . 

For A > 0, pinned interfaces at = Fc are self-affine 
with a roughness exponent a ~ 0.63. Moving interfaces 
are not self-affine which we demonstrated by measuring 
the correlation functions C'^{r,t) and c(r) (see Figs. 2- 
4). This is in comparison to the simpler behavior of Eq. 
(^) with A = 0, where not only pinned interfaces are self- 
affine (on length scales I <C ^) but also moving interfaces 
on time scales r ^ The latter can be seen from the 
fact that the height-height correlation function c(r) for 
an interface moving with constant velocity scales with 
the same growth exponent (3 ~ 0.88 as the global width 
at f = Fc (see Fig. 3a and Table I). The values for the 
exponent /?,„ (Eq. (|^)) are obtained by the scahng of the 
global width w{t) on time scales t and correspond 

to those of the Edwards- Wilkinson equation (A = 0) ]3^ ] 
and the KPZ equation (A > 0) ||l^, respectively. 

The results for the critical exponents characterizing the 
depinning transition of Eq. (^ with A > are summa- 
rized in Table II. 



Exponent 


DPD 


continuum eq. 


automaton 


roughness, a 
growth, P 
correlation length, v 
velocity, 9 


0.633 
0.633 
1.733 
0.636 


0.63 ± 0.01 
0.64 ± 0.02 

1.7 ± 0.1 
0.64 ± 0.05 


0.63 ± 0.01 
0.64 ± 0.02 
1.72 ± 0.03 
0.63 ± 0.02 



Table II. Comparison of the critical exponents of the 
DPD model with our numerical results. 

The numerical results for both, the continuum equa- 
tion (fol) and the automaton are in excellent agreement 
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with the DPD modeL We therefore expect that the criti- 
cal exponents of the anisotropic depinning model Eq. (|p) 
are exactly given by the exponents of directed percola- 
tion. 

In the DPD model, the dynamic and the static behav- 
ior is determined by directed percolation paths of pinning 
sites. Due to the restricted solid-on-solid condition of the 
directed percolation paths, the pinned regions of the in- 
terface have small gradients (bounded by a slope one) 
|p7| . In contrast, the laterally moving regions have a lin- 
ear slope of about two This behavior is analogous 
to that of Eq. regions of the interface with large 

slopes are likely to move due to the positive contribution 
of the term A(V/i)^ to the driving force. Regions of the 
interface with small gradients, the other hand, are easier 
to pin, due to the smaller contribution of A(V/i)^. 

The observation that the interface motion is mainly 
due to the gradient term A(V/i)^, causes a clustering of 
growth sites, which can be understood as follows. A mo- 
tion of an interface element h{x) h{x) + dh increases 
the contribution of A(V/i)^ to the local force at h{x + dx) 
or h(x — dx). Thus, this neighboring interface element is 
likely to be the next new growth site, resulting in a cluster 
of growth sites to be formed. The moving regions have 
larger slopes than the pinned parts. As a consequence, 
a moving interface is not self-affine. The clustering of 
growth sites in Eq. ^ does not destroy the dynamical 
scaling of global quantities of the interface, such as the 
width and the velocity (see Eqs. (||-§|)). The scaling re- 
lation (||) is also fulfilled by the exponents of the DPD 
model and by our numerical results. 

The behavior of moving interfaces deserves further in- 
vestigation to understand better the concepts of scaling 
and self-affinity (see Sec. IV B). Further, it would be very 
interesting to construct a rigorous proof that Eq. (||) and 
the DPD model are in the same universality class. 
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